l/f a fluctuations in a ricepile model 
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The temporal fluctuation of the average slope of a ricepile model is investigated. It is found that 
the power spectrum S(f) scales as 1/ f a with a ~ 1.3 when grains of rice are added only to one end 
of the pile. If grains are randomly added to the pile, the power spectrum exhibits l// 2 behavior. 
The profile fluctuations of the pile under different driving mechanisms are also discussed. 
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The term flick noise refers to the phenomenon that a 
signal s(t) fluctuates with a power spectrum S(f) oc 1/ f a 
at low frequency. Since the exponent a is often close to 
1, flick noise is also called l/f noise. The power spec- 
trum of a signal is defined as the Fourier transform of 
its auto-correlation function C(to,t) = {s(to)s(t)), where 
(■ • •} denotes ensemble average. Usually the signal s(t) 
under study is stationary and its auto-correlation func- 
tion only depends on t— to- The auto-correlation function 
can be alternatively calculated as 



C(r) 



?(t)s(t + r)dt. 



(1) 



If the signal is real-valued, its power spectrum is just 

S(f) = \s(f)\ 2 (2) 

according to Wiener-Khinchin Theorem. Here s"(/) is the 
Fourier transform of the signal 



s(f) 



s (ty 2 * ft dt. 



(3) 



In nature and laboratory, many physical systems show 
flick noise. For example, flick noise appears in a variety 
of systems ranging from the light of quasars [Q to water 
flows in rivers [01, music and speech B, and electrical 
measurement ||||. Despite its ubiquity, a universal ex- 
planation of the flick noise is still lacking. Bak, Tang, and 
Wiesenfeld have proposed that self-organized criticality 
(SOC) may be the mechanism underlying flick noise. 
They also demonstrated their idea of SOC with a cellular 
automaton model, the BTW sandpile model. However, 
the model they proposed does not exhibit flick noise in its 
ID and 2D versions f|. Especially the ID BTW sandpile 
model does not even exhibit SOC behavior. Experiments 
|| on piles of rice grains were done to investigate if real 
granular systems display SOC behaviors. It was found 
that ricepiles with elongated rice grains exhibit SOC be- 
haviors. The key ingredient here is the heterogeneity 
of the system. Because elongated rice grains could be 
packed in different ways ||, the ricepile thus has many 
different metastable states. A ID model of ricepile was 
proposed by the same group of authors; avalanche dis- 
tribution and transit-time distributions were discussed 



||. The model was later investigated by several other 
authors including us ]To|-^2]|. Attentions were paid to 
the avalanche dynamics and transit-time statistics of the 
model. In this paper, we are concerned with temporal 
fluctuation of the average slope of the system. We want 
to investigate if it has flick noise or flick-like behavior. 

The ricepile model is defined as follows: Consider a 
one-dimensional array of sites 1, 2, • • ■, L. Each site con- 
tains an integer number hi of rice grains. Here hi is called 
the local height of the ricepile at site i. The system is 
initialized by setting hi — for all sites. This means the 
ricepile is to be built up from scratch. The system is 
driven by dropping rice grains onto the pile. If one grain 
of rice is dropped at site i, then the height of column i 
will increase by 1, hi — ► hi + 1. With the dropping of rice 
grains, a ricepile is built up. The local slope of the pile 
is defined as Zi = hi — /ij+i. Whenever the local slope Zi 
exceeds certain threshold zf (specified in the following), 
site i will topple and one grain of rice will be transferred 
to its neighbor site on the right. That is, hi — > hi — 1 and 
hi+i — * hi+i + 1- In this model, rice grains are allowed 
to leave the pile only from the right end of the pile, while 
the left end of the pile is closed. So the boundary con- 
dition is kept as ho — hi and hr,+i = 0. The ricepile is 
said to be stable if no local slope exceeds threshold value, 
that is, Zi < z\ for all i. Rice grains are dropped to the 
pile only when the pile is stable. In an unstable state, all 
sites i with z\ > zf topple in parallel. The topplings of 
one or more sites is called an avalanche event, and the 
size of an avalanche is defined as the number of topplings 
involved in the avalanche event. 

The values of zf s are essential to the definition of the 
model. As in Ref. jnj, every time site i topples, the 
threshold slope zf will take a new value randomly chosen 
from 1 to r with equal probability. Here r is an inte- 
ger no less than 1. The parameter r reflects the level of 
medium disorder. The larger the r, the higher the level 
of medium disorder of the system. If r = 1 the model 
becomes the ID BTW sandpile model. When r = 2, the 
model reduces to the model studied in reference In 
ref. 1 10 1, we have investigated the effect of disorder on the 
universality of the avalanche size distribution and transit 
time distribution. 
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In the present work, we performed extensive numerical 
simulations on the system evolution. Let us first study 
the case where rice grains are added only to the left end 
of the pile, i.e., only to the site i = 1, which is in accor- 
dance with the experiment setup ||. We shall refer to 
this driving mechanism as the fixed driving. When the 
stationary state is reached, we record the average slope 
z(t) = hi(t)/L of the pile after every avalanche. Here 
time t is measured in the number of grained added to 
the system. Thus we got a time series z(t). Typical re- 
sults about the fluctuation of the slope z(t) can be seen 
in FigJ^.a. We calculate the power spectrum of this time 
series according to Eq. (||) . We find that for not too small 
systems the slope fluctuation displays l//-like behavior. 
In fact, we got a power spectrum S(f) scaling as 1/ f a at 
low frequency, with the exponent a = 1.3 ±0.1. 

We note the following points: 

(a) the temporal behavior of the ricepile model is dra- 
matically different from the ID BTW sandpile model. In 
the ID BTW model, the system has a stable state [||, 
which, once reached, cannot be altered by subsequent 
droppings of sand grains. Thus the ID BTW model does 
not display SOC behaviors. The slope of the ID sandpile 
model assumes constant values. The power spectrum is 
thus of the form S(f) — 5(f). In words, there is only 
dc component in the power spectrum for the ID BTW 
model. 

(b) With the introduction of varying threshold zf into 
the ricepile model the behavior of the system becomes 
much more rich. The avalanche distribution follows 
power law Jl0[ , which is an important signature of SOC. 
Besides that, the temporal fluctuation of z(t) has a power 
spectrum of flick type at low frequency. 

We notice that the power spectrum is flattened at low 
frequency for small system. This is a size effect. It is be- 
lieved that flick noise fluctuation is closely related to the 
long range spatial correlation in the system. For small 
system, the spatial correlation that can be built up is 
limited by the system size, thus the long range temporal 
correlation required by flick noise is truncated off at low 
frequency. Our numerical results verified the above dis- 
cussion. In Fig. ||, we show the power spectrum of z(t) 
for different system sizes. It is clear that when the sys- 
tem size increases the 1//" behaviors extends to lower 
and lower frequency. 

We also investigated the effect of the level of disorder 
on the power spectrum. This was done by simulating the 
system with different values of r. In Fig. || we show the 
power spectra for the system with different r. It can be 
seen from this figure, that the power-law (straight in the 
log-log plot) part extends to lower frequency for higher 
value of r. The effect can be understood by the following 
discussion. When the level of disorder is increased, the 
amplitude of the slope fluctuation also increases. Larger 
amplitude fluctuation of z(t) requires more grains to be 
added to the system. Hence the period of this fluctu- 
ation increase, which gives rise to the increase of low 
frequency components in the power spectrum. We have 



made statistics on deviation of the ricepile slope from 
its average value in the stationary state. We found that 
the greater the r, the greater the deviation. Note that 
the value of r > 1 does not affect the exponent a, as it 
does not alter the universal avalanche exponent r for the 
avalanche distribution |l(J . 

It is also interesting to study the effect of driving mech- 
anism on the temporal behavior of the system. Now, in- 
stead of dropping rice grains to the left end of the pile, we 
drop rice grains to randomly chosen sites of the system. 
This way of dropping rice grains represents a different 
driving mechanism, which we shall refer to as the ran- 
dom driving. Numerical simulations with this random 
driving were made and the time series z(t) was recorded. 
For this way of driving, typical result about the fluctua- 
tion of z(t) can be seen in Fig. |l|.b. It seems that z(t) 
now fluctuates more regularly than it does for the fixed 
driving. We found that under random driving the tem- 
poral fluctuation has a trivial l// 2 behavior as in the 
case of various previously studied sandpiles . In Fig. 
|3|, we compare the power spectra of the two cases with 
different driving mechanisms. Two groups of curves are 
shown in this figure, each with a different exponent a. 
For the fixed driving, we have a ~ 1.25, while for the 
random driving we have a w 2.0. The distinction of the 
two behaviors is quite clear. This shows that the driving 
mechanism, as an integrated part of the model, has an 
important role in the temporal behaviors of the system. 
Recent work on the continuous version of BTW sandpile 
model HH, and that on quasi- ID BTW model jSj also 
showed that certain driving mechanisms are necessary 
conditions for these models to display 1// fluctuation for 
the total amount of sand (or say, energy). 

Here we present an heuristic discussion on why the ex- 
ponent a is smaller for fixed driving than for random 
driving. For fixed driving, the height hi at site i = 1 
changes almost every dropping, and so is the average 
slope z(t). Thus the high frequency component of the 
z(t) fluctuation has a heavier weight in its power spec- 
trum, and this will make the exponent a smaller. While 
for the random driving, the rice grains are dropped to the 
pile at random sites, the spatial correlations previously 
built up can be easily destroyed, making the z(t) behave 
more or less as an random walk. So the exponent a for 
this case shall be very close to 2. For random driving, 
each site has the same probability in receiving a grain in 
every drop. For the fixed driving, however, only the left 
end site receive grains, so there is in some sense breaking 
of symmetry, which would lead to different scalings in the 
slope fluctuations. 

To see more about the effect of driving mechanism, it 
is helpful to investigate the profile fluctuations of the pile 
under different drivings. Because avalanches change the 
surface of the ricepile, the pile fluctuates around an av- 
erage profile, and the size of the fluctuations characterize 
the active zone of the ricepile surface. Let the standard 
deviation of height at site i be (Jh(h L) = \J (tif) — (hi) 2 . 
Here (• • •} represents average over time. Then we calcu- 
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late the profile width of the ricepile, w = \ YN (7h{i, L), 
which is a function of the system size. In Fig. |4|, we show 
the profile width of the ricepile for different r and differ- 
ent driving mechanisms. It can be seen that w scales with 
L as w (x L x , and that x = 0-25 ± 0.01 for fixed driving, 
X = 0.09 ± 0.01 for random driving. For given driving 
mechanism, w increases with increasing r. As we stated 
in Ref . jl(j , the parameter r reflects the level of medium 
disorder in the rice pile. Although the level of disorder 
does not change the scaling exponent x, it does affect 
the amplitude of fluctuations. For greater r, the profile 
width is larger. It can be seen in the figure, that the data 
points for r = 4 are above that for r = 2, for given driv- 
ing mechanism. A recent experiment showed that piles 
of polished rice gra ins has a smaller profile width than 
unpolished ones pa], which have higher level of medium 
disorder. 

In conclusion, we have studied the temporal fluctua- 
tion of the slope of ricepile model in its critical stationary 
state. The power spectrum of this fluctuation is closely 
related to the driving mechanism. When the rice grains 
are dropped to the left end of the pile, the slope fluctu- 
ates with a flick-type power spectrum, with the exponent 
q = 1.3 ± 0.1. When the driving mechanism is changed 
to random driving, the model displays l// 2 behaviors. 
Greater system size and higher level of disorder will ex- 
tend the frequency range where the power spectrum has 
the form 1// Q . 
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FIG. 1. The time evolution of the average slope z(t) of 
the ricepile after the stationary state has been reached. The 
parameters used for this figure are L = 20, r = 4. (a) rice 
grains are added to the pile at site i = 1. (b) rice grains are 
added randomly to the pile. 
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FIG. 2. The power spectrum of z(t) for different system 
sizes, with fixed driving and r = 2. The curves (from top to 
bottom) are for L — 40, 60, and 100 respectively. The dashed 
line is a curve for y oc x~ 1,25 for reference. The number of data 
points for Fourier transformation is N = 2048, and the results 
are obtained by averaging over 500 samples. f c = 1/(2A) is 
the Nyquist critical frequency, where A is the time interval 
between two successive points in the Fourier transformation. 
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FIG. 3. The power spectrum of the slope fluctuation for 
a ricepile with L — 100. Two groups of curves are shown 
in this figure, with different exponents a. For reference, the 
upper dashed line is a curve for y oc x~ 1,25 , while the lower 
dotted line shows a curve y oc x~ 2 . In each group of curves, 
the upper curve is for the case r = 4, and the lower one is for 
r = 2. 
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FIG. 4. The profile width of the ricepile scales with the 
system size. Data from numerical simulations are shown in 
open (for r = 2) or filled ( for r = 4) symbols. Solid lines 
are power law fit to the data for fixed driving, with power 
law exponent \ — 0-25 ± 0.01, and the dashed lines are for 
random driving, with x = 0.09±0.01. For every run of numer- 
ical simulation, statistics were made over at least 10 s grain 
dropping. 
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